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Abstract 

Full radiative transfer in 3D space has been solved to pursue the reion- 
ization history in an inhomogeneous universe. It has been shown that the reion- 
ization of an inhomogeneous universe is not a prompt event, but a fairly slow 
process. Also, QSO absorption hne systems are simulated with using the results 
of radiative transfer calculations. Lya continuum depression implies that the 
metagalactic UV intensity decreases rapidly with z at z > 5, and the reionization 
must have taken place between z = Q and 10. Finally, it is stressed that Ha 
forest is a more powerful tool to probe the reionization history and the density 
fluctuations in the universe at 2; > 5. 

1. Introduction 

The cosmic reionization is one of the most significant issues in cosmology, 
which is closely related to the formation of QSOs and galaxies. The information 
on the ionization states in the universe has been accumulated by the observations 
of Lya absorption lines in high redshift QSO or galaxy spectra. 3D cosmological 
hydrodynamic simulations (Cen et al. 1994; Miralda-Escude et al. 1996; Gnedin 
& Ostriker 1996; Zhang et al. 1997) have revealed that the Lya absorption sys- 
tems can be accounted for in terms of the absorption by intergalactic density 
fluctuations. However, all of these works have been based upon optically-thin 
or local optical-depth approximations. Recently, it has been shown that the ra- 
diative transfer effects of ionizing radiation could strongly affect the ionization 
structure (Razoumov & Scott 1999; Gnedin 1999; Nakamoto et al. 1999). Madau 
(1995) and Haardt & Madau (1996) as well have considered cosmological radia- 
tive transfer on the assumption of ID semi-inflnite slab. Here, we present 3D 
radiative transfer calculations on the cosmic reionization and the formation of 
QSO absorption line systems. 
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2. Method 

To reproduce an inhomogeneous universe, we generate density fluctuations 
based upon the Zel'dovich approximation. In the present calculations, we assume 
a standard cold dark matter cosmology, i.e., ^cdm = 0.95, OBaryon = 0.05, erg = 
0.6 with the Hubble constant of 50 km s^^ Mpc~-^. 

We calculate the ionization degree by assuming the temperature of lO'^K 
and ionization equilibrium. The simulation box is irradiated by the isotropic 
UV background radiation of a power law- type spectrum, /q = /2ilO~^^ erg cm~^ 
s^^ Hz~^ sr~^. The so-called "proximity effect" of Lyman alpha absorption hnes 
requires the diffuse UV radiation to be at a level of /21 ~ 1 at 2 < z < 4 (Bajthk, 
Duncan, & Ostriker 1988; Giallongo ct al. 1996). 

The UV radiation fields are obtained by solving the three-dimensional 
steady radiative transfer equation. For the implementation of the radiative trans- 
fer on a massively parallel supercomputer, we have developed a new scheme 'Se- 
quential Wave Front Method'. We have used 128^ grids in space, 128^ in direc- 
tions, and 6 for important line frequencies of hydrogen and helium. Since the con- 
tinuum is analytically integrated, the calculations are quite accurate even if just 
6 frequencies are solved. The total operations amounts to about 1 Tflops*hour. 
The calculations have been performed on the CP-PACS in University of Tsukuba. 

3. Results &6 Discussion 

3.1. Cosmic Reionization 

Around z — 15, underdense regions are reionized and percolate, leaving 
overdense neutral islands due to the self-shielding effects. Since the absolute den- 
sity decreases with time in a linear regime of density fluctuations, the reionization 
proceeds and the ionized sea encroaches onto the neutral islands, leaving fllamen- 
tary self-shielded regions at 2; ~ 9. The highly ionized regions expand further 
at z ~ 7, but the universe is not perfectly transparent against UV radiation in 
the sense that the collective optical depth considerably reduces the incident ra- 
diation. Around z = 5, the universe becomes transparent against background 
UV and the overall reionization has been fulflUed. The simulations have revealed 
that, owing to the radiative transfer effects, the reionization in an inhomogeneous 
universe is not a prompt event , but fairly slow process, in contrast to the prompt 
reionization in a homogeneous universe. 



3 



3.2. Lya Absorption Lines and Continuum Depression 

The resultant ionization degrees in the universe are different from place 
to place by more than three orders of magnitude. Due to such inhomogeneous 
ionization structure, relatively low ionization regions could produce strong ab- 
sorption in quasar spectra. To make a direct comparison with the observations, 
we simulate absorption lines. First, we focus on Lya absorption. To match the 
recent observations by the Keck telescope, we adopt the resolution of R=45000 
and the variance of 0.04, and assume the Voigt profile of lines. The simulated 
Lya absorption features are shown in Figure 1. 



Fig.l - The simulated Lya absorption lines against 
wavelength at z = 3 [top panel) with I21 = 1, z = 
4 (second) with I21 = 0.1, and z = 5 (third) with 
721 = 0.01. The bottom panel is the diagram of Lya 
continuum depression (thick gray curves) against red- 
shifts. Symbols arc observations. Also, the mean neu- 
tral fractions Xfjj (thin curves) are shown. The same 
line types correspond to the same UV intensity. 



Fig.2 — Same as Fig. 1, but for the Ha absoiirtiori. 
There are no observation data so far for Ha absorption 
at high redshifts 



In order to compare quantitatively the simulations with observations, we 
have assessed the so-called continuum depression. Da- Figure 1 shows that any 
model with a constant UV intensity does not match the observed trend that 
Da tends to grow quickly at higher redshifts up to 5. This implies that the 
metagalactic UV intensity must decrease rapidly with z at z > 5 hy two orders 
of magnitude at least. If the well fitted value of I21 = 0.01 at 2 = 5 is unchanged 
also at higher redshifts, the reionization epoch is estimated tohe z ^ 9. If the UV 
intensity decreases in a similar fashion also at higher redshifts, the reionization 
epoch is 2; ~ 6. Thus, it is concluded that the cosmic reionization must have taken 
place between z — 6 and 10. However, the Lya absorption is not appropriate 
to determine the reionization epoch more accurately, because, as seen in the 
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absorption features in Figure 1, Lya is too strongly depleted even if the mean 
neutral fraction is less than 10~^. In other words, Da is no longer sensitive to 
Xhi around the cosmic reionization epoch. 



3.3. Ha Forest 

As shown in the previous subsection, Lya has too high line opacity to 
probe the universe at z > 5. Therefore, three conditions are required for a line 
in order to investigate the universe at 2 > 5: (1) it has lower line opacity than 
Lya, (2) line emission is detectable, and (3) it has lower extinction against dust 
because young star-forming galaxies are often dust-enshrouded. The most favor- 
able solution is Ha absorption lines. In Figure 2, the Ha absorption features and 
continuum depression are shown. The Ha absorption is relatively weak around 
2 = 3, while it is sensitive to the ionization degrees at 2 > 4. Also the Ha 
continuum depression traces the reionization history more accurately. Therefore, 
it is concluded that Ha forest is a more powerful tool to probe the universe at 
2; > 5. Ha forest has been never detected so far. The reason comes from the 
fact that Ha has much weaker opacity than Lya. From observational points of 
view, the continuum depression can be detected by low- dispersion spectroscopy 
or narrow-band photometry. Furthermore, Ha forest is subject to less UV bump 
effects for AGNs compared to Lya forest. The wavelengths of Ha forest drop on 
3^m < Ana < 7/im at 4 < 2; < 10. Thus, the observations can be done with Subaru 
IRCS, IRIS, SIRTF, NGST, or H2/L2. If one can obtain the absorption features 
with the resolution greater than 10000, one can recover the density fluctuations 
at high redshifts. They allow us to determine the linear amplitude of pregalac- 
tic perturbations which is by no means measured in the GBR due to the strong 
Sunyaev-Zeldovich effects in galactic scales. If one has the amplitude of linear 
density fluctuations at galactic scales, one can not only set the initial condition 
for galaxy formation, but also make more reliable determination of cosmological 
parameters. 
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